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Currently, spin tunneling at very low temperatures is assumed to proceed as an incoherent 
sequence of events that take place whenever a bias field h(t) that varies randomly with time t 
becomes sufficiently small, as in Landau-Zener transitions. We study the behavior of a suitably 
defined coherence time t$. Coherence effects become significant when > Th, where th is the 
correlation time for h(t). The theory of tunneling of Prokof'ev and Stamp (PS), which rests on the 
assumption that < th, is extended beyond this constraint. It is shown, both analytically and 
numerically, that > th when ThSh < h, where 5h is the rms deviation of h. Equations that give 
t,/, and the tunneling rate as a function of ThSh both for ThSh > h, where the theory of PS hold, and 
for ThSh < h, where it does not, are derived. 
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Magnetic quantum tunneling (MQT) of large spins has 
become a subject of much interest. Much of it stems from 
the expectation that MQT will contribute to our under- 
standing of decoherence. El Decoherence effects play an im- 
portanLrole in the transition from quantum to classical 
physicso and in quantum computing!] Magnetic quan- 
tum tunneling (MQT) through thermally activated states 
is thought to proceed incoherentlyoEl However, at suffi- 
ciently low temperatures,|4:unneling must proceed though 
the ground state doublets. Two-state systems aig-often 
used to model tunneling under such conditions.Q'El Ex- 
isting theory treats tunneling between these two states 
as an incoherent sequence of attempts that take place 
whenever the difference between the two-state energies 
becomes sufficiently smalLta It is in the spirit of Landau- 
Zener schemes that are often used in the theory of non- 
adiabatic transitions to estimate transfer rates between 
two states whose energies sometimes cross.E2l 

Here, I study coherence effects in MQT by means of 
the following two-state model, 



H = ea z + Acr x — h(t)a z 



(1) 



where Ti. is the Hamiltonian, a z and a x are Pauli spin 
matrices, e is an energy that is constant in time, h(t) 
varies with time t following a stationary random Gaus- 
sian process, and A stands for the tunneling frequency 
(H = 1, and k = 1, where k is Boltzmann's constant, 
from here op.)_ Two-state systems have been used to 
study MQTj!r@ but they have, also been used to study 
decoherence in other .systems, til such as in (hypothetical) 
quantum computers.0 The heat bath is replaced by h{t) 
in this model. Consequently, the model is applicable only 
if A <C T and e«T, where T is the temperature. This 
is amply satisfied in present day experiments in MQT 
when no transverse field is applied. 

The aim of this paper is to study within the model de- 
fined above the behavior of the coherence time 70 (defined 



below) and of the tunneling rate T as a function of ThSh, 
where Sh is the rms deviation of h(t) and Th is its corre- 
lation time. Unless stated otherwise, t/jA <C 1, A< 5h, 
and A <C (ThSh)Sh are assumed to hold everywhere. It 
is shown that decoherence effects are partially averaged 
out when ThSh < 1. I derive, and support with numer- 
ical results, that if ThSh < 1, then t^/th ~ l/(ThSh) 2 , 
instead of ps 1/Sh that is shown to hold if ThSh > 1. 
Furthermore, instead of the relation. 
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derived below for the tunneling rate under conditions 
of full incoherence (that is, for < r^), which is the 
underlying assumption in the spin tunneling theory of 
Prokof'ev and Stamp,El the relation 
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where f, = 2t> 1 (<5/i) 2 , is derived under the condition 
that ThSh < 1. Note that « t$ if ThSh < 1. Finally, 
concluding remarks about the likely hood that ThSh < 1 
is realized in Fes & n d Mni2 are made. j—. 

In order to define a measure of coherence JL3 consider 
the probability p a {t) that the system evolve in time t 
from an intial state into a final state |er). We can 
then write, 

P „(t+T)= Y,™&°',<r"-,T)wm)){mw"), (4) 



where w^ctW) = (xr"\U^{{h},r)\a){a\U{{h},T)\a'), 
U({h},T) = T expH Jo H(tx)dti] for a given h(t) in Eq. 
(1), and T is a time ordering operator. Averaging Eq. 
(i) over different h(t) time sequences, then when Th for 
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h(t) fulfills Th <C r, t, its right hand side breaks up 
into sums of products of two terms, (w(a,a' ,er";r))/i 
and ((a'\ip(t))(i(!(t)\a"))h, where (}/, stands for average 
over all different histories of h(t). We are interested 
here in establishing the physical conditions under which 
(w(a,a' ,a";r))h — > £ -/ >ct /'W ct ,ct'( t ) ensues, because this 
leads to 



2><r(i + T) =^W CTjff '(r)^(t). 



(5) 



No interference terms from Eq. (|j) appear in this equa- 
tion. Thus, probabilities, rather than probability ampli- 
tudes, for alternative paths to a final state are summed. 
Its applicability is,therefore, tantamount to incoherence. 
Solutions to Eq. Q5J) exhibit no oscillations. For any two 
state system, they relax to the final state exponentially. 
Consider the off-diagonal quantity, 



(w(a, a',a";T)} h 



(6) 



and 



for a' ± cr", where Pl = (\{a"\W {{h}, r)\a) 
Vi = (IWJU ({h} , t)\ct')\ 2 ) h- As shown above, q — leads 



to Eq. (|5j), and therefore to incoherence. On the other 
hand, the system evolves in time in a pure state if there is 
no averaging over different h(t) time evolutions. No de- 
coherence occurs then, and it follows from the definition 
of q that | q \= 1 then. Quantity q provides, therefore, 
an appropriate measure of coherence between the two 
alternative paths a' and a" . This is in analogy to an in- 
terference experiment that may be performed to test the 
degree of coherence between the parts of a particle-wave 
that go through two slits placed at a' and a" . Let the co- 
herence time be defined as the smallest time such that, 
say, | q(r) |< 1/2 for all r > ta,. (Below, q is independent 
of cr, a' , a" throughout. Accordingly, these variables of q 
are omitted from here on.) 

Numerical results are obtained for the system's time 
evolution for a given h{t) time sequence by iterating a 
finite difference Schrodinger's equation, 

Ht + Dt) = {1 - iDtU{t) - [DtH(t)] 2 /2}iJj(t), (7) 

where Tt is defind in Eq. (Q), ip is a two-component 
state function, and Dt is a suitably small time interval 
(Dt < 10 _6 /A keeps the normalization condition satis- 
fied within 1 part out of 1000 throughout the reported 
simulations). 

The values of h(t) are assumed to follow from a sta- 
tionary random Gaussian process, with a distribution for 
h centered on 0, and h(t) is correlated over time ac- 
cording to (h(t)h(0))t — (Sh) 2 exp(— t/r h ). This is ac- 
complished with Langevin's (finite difference) equation, 
h(t + Dt) = h(t) - Dt[h(t)/T h + f(t)}, where /(*) is as- 
signed an independent random number at each value of 
t. All values of r h and Dt used fulfill Dt < r h /10. The 
values of h obtained follow the prescribed behavior for it 
if (/) = and (/(t)/(t')) = 5 t , t ,(2-Dt/T h )(Sh) 2 /(T h Dt), 



where S^t' is the Kronecker delta function. The correc- 
tion term Dt/r^ comes about from applying the fluctua- 
tion dissipation theorem to the finite difference Langevin 
equation used here. 

Numerical results obtained for the tunneling probabil- 
ity and for q are shown in Fig. 1 for Sh/A = 50 and 
ThA = 10~ 4 . Note that even though 5h is 50 times 
larger than A, the tunneling probability undergoes os- 
cillations, and, correspondingly, 1 only after a time 
that is much larger than A -1 . As shown below, the 
high degree of coherence exhibited in Fig. 1 occurs when 
Th(Sh) 2 < A. It may be realized in MQT experiments 
by enlarging A, .through the application of an external 
field, as in Refs.,E3o where A « 30K is several orders of 
magnitude larger than the linewidth [see Eq. below] . 




FIG. 1. The probability p that tunneling occurs in time t 
and the off-diagonal quantity q are shown versus tA. The 
data points shown follow from an average over 4 x 10 3 time 
evolutions with different realizations of h(t), with a vanishing 
time average value of h, Sh/A — 50, and r^A = 10 -4 . 



This extreme condition is worth considering, because it 
serves to illustrate how decoherence effects can be drasti- 
cally supressed. Quite simply, rapid fluctuations in h are 
almost completely averaged out. More specifically, they 
are scaled down by the factor \/ /r^ in time r^. Thus, 
5h/A — 50 is scaled down by 100 in Fig. 1. 

Numerical results that illustrate the effects of partially 
averaging decoherence out are exhibited in Fig. 2. 
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FIG. 2. Tunneling probability p versus time. The curves 
shown follow from an average over 4 x 10 3 time evolutions 
with different realizations of hit), e = 0, Sh/A = 10 3 , and 
thA = 10~ 5 . From data points obtained for q (not shown), 
the value of time where q = 1/2, that is t^, can be obtained. 
It turns out that ~ 0.07/A. The full line is for p; the dotted 
line is for the prediction that follows for p assuming incoher- 
ence for times shorter than tk\ the dashed-dotted line is for 
(At) 2 , that is fully coherent evolution of p; finally, the long 
dashed line is for the assymtotic evolution after decoherence 
sets in. 

To start the derivation of the tunneling rates, 
consider the solution to Schrodinger's equation, 
yj(t) = U ({h},t)Ui({h},t)il>(0), where U ({h},t) = 
cxp[-i J Q Ho(T)dT], Ho(t) = [e- h(r)]a z , Ui({h},t) = 

Texp [-*/„* Vj(T)dT], Vj(t) = ul({h},t)VU {{h},t), 
and V = Ao~ x . The probability amplitude a(t) that 
tunneling take place in time t, that is, that a transition 
from state a z — — 1 to state o z = +1 occurs, is 

| a(t) |=| f dt 1 (llAo-xe^/o^oWirj _ x) ^ (8) 



in first order perturbation theory, where use has been 
made of the relation exp(Aa z ) a x — o~ x exp(— Aa z ) for 
any c-number A. Now, since H.o(t) | ±1) = ±[e — h(r)] \ 
±1), it follows that g l H Q {T)dT can be replaced in Eq. 
(|) by [e - where H{t) is defined by tH(t) = 

r Q dr'h{r'). 

In order to get some feeling for the physics behind 
the results to be derived, a rough argument that leads 
to semiquantitative relations for the coherence time and 
the tunneling rate V is given first. For shortness, the 
argument is restricted to the condition e = 0, but the 
derivation itself of Eqs. (0) and (^) is not. It follows 
from Eq. (^), and the definitions of q and of that, in 
general, 



r « 

where 5H(t$) is the rms value of H(t$). 
SH( Trj) ) pb Sh if r h Sh > 1. Therefore, 



(9) 



if ThSh > 1. There is then complete incoherence, pthat 
i s > T 4> T h- The scheme of Prokof'ev and StampO, of 
an incoherent sum of tunneling probability increments 
that occur whenever h(t) becomes comparable to A, is 
then justified. If, on the other hand, ThSh < 1, then 
5H(t^) ss Sh/ y/r^/Th, since t^/t/j is, in a rough sense, 
the number of times that h fluctuates independently dur- 
ing time 70. Therefore, substitution into Eq. (g) gives 



t^/th sa \/(T h 8Ky 



(11) 



if Tft(5/i < 1. This equation exhibits explicitely how av- 
eraging out of fluctuations makes larger than if 
77i<5/i < 1. Data points obtained numerically for are 
shown in Fig. 3. Note that Eq. ( pd| ) also implies that co- 
herence times become larger than tunneling times when 
T h 5h 2 < A. 
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FIG. 3. Data points from numerical time evolutions with 
e = for t^/th and for for TSh/A versus ThSh. ■ and ♦ 
stand for Sh = 1000, o, •, and □ stand for 300, 75, and 30, re- 
spectively. Each data point for t^/th and for TSh/A, stands 
for an average over at least 10 3 and 10 4 h(t) time evolutions, 
respectively. The straight lines for t^/th are guides to the 
eye, with slopes 1 and 2; the straight lines for TSh/A follow 
from Eqs. (g|) and (^). All data points shown are for values 
of r h that fulfill r h A < 0.1 and T h (Sh) 2 / A > 1. Error bars 
are covered by the symbols for the data points. 

The tunneling rate T is estimated next for e = and 
T hSh < 1. It follows from Eq. (^) that the probabil- 
ity p[t, H(t)\ for a given H (t)] that a spin has tunnelled 
at time t is p[t,H(t)} w A 2 i 2 for tH[t) < 1. It may 
seem strange that a spin under a large field tunnels as 
fast as one for which H = even though the line shape 
is known to be a Lorcntzian function, but note that 
A 2 (A 2 + H 2 )' 1 sin 2 (fit), where ft 2 = H 2 + A 2 , is ap- 
proximately A 2 t 2 for fit < 1, wich is independent of H. 

2^2 



Furthermore, Averaging over all h(t) sequences leads to p(r^) « A 2 t^ 
Equation (g) gives p(t) f=s p(r^)t/T^ for < > t^, since 
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T h % Tfj, under the assumed condition ThSh < 1. It fol- 
lows then that, T ~ A 2 /[T h {Sh) 2 ] for s = 0. 

The derivation of Eqs. (|), ©, ©, and @, from 
which we digressed below Eq. (^) is now resumed. It is 
assumed that h(t) follows a stationary random Gaussian 
process, centered on [a non-zero value simply redefines e 
in Eq. (1)], and h(t) is correlated over time, with a corre- 
lation time Th defined by Th{h 2 ) = J Q (h(t')h(t' + t))hdt, 
where (. . )h stands for an average over all h(t) histories. 
The probability, 

P(t) = A 2 [ dt 2 f 2 dh{e 2l ^ [h{T) - £]dT + cc}, (12) 
Jo Jo 

follows straightforwardly from Eq. (||) and the statement 
following it. Now, J^ 2 h(r)dT follows a Gaussian distri- 
bution since h(t) follows, by assumption, from a station- 
ary random Gaussian process. Therefore, after averaging 
over all histories of h(t), Eq. ( |l2] ) becomes, 

P(t) = A 2 f dt 2 [ dt x e - 2F{tl '^ cos [2e(t a - h)}, 
Jo Jo 

(13) 

where F(t\,t 2 ) = dr x dT 2 (h(T 1 )h(T 2 ))h- Now, 
F(h,t 2 ) ~ 2{8h) 2 T h \t 2 - ii| if r h < \t 2 - and a bit 
of reflection shows that this relation applies if rSh <C 1 
and t 3> \/[Th{8h) 2 ] (that is, if t 3> ta), which leads to 
Eq. (§j). On the other hand, F(h,t 2 ) ~ (<5/i) 2 |i 2 - *i | 2 if 
Th ^ |*2 — ii|, and this relation is applicable if ThSh » 1 
and t 3> \/{Th5h) (that is, if t ^> r^), which leads to Eq. 
(|J). Similarly, 

<- ± jmf.** m '''™ M ' (14) 

follows from Eq. (^|) and the assumption of a stationary 
random Gaussian h(t). By the same arguments that led 
to Eqs. (Q), and (^) above, equations ([!(]) and ( pd] ) follow 
straightforwardly, independently of e. 

Numerical results obtained for T and for are shown 
in Fig. 3. There is good agreement with Eqs. 

In conclusion, coherence effects have been shown to set 
in when the coherence time is longer than the correla- 
tion time Th for h(t), and that this comes about because 
decoherence effects that arise from h(t) partially average 
out when Th < 1/Sh. Coherence times are then larger 
than Th- Equations (^|) and (|§|) which give and the 
tunneling rate as a function of ThSh both-.for ThSh > H, 
where the theory of Prokof 'ev and StampH hold, and for 
ThSh < h, where their theory does not hold, are derived. 

The condition ThSh < 1, which gives > Th, would be 
fulfilled in spin tunneling systems in which nuclear spin 
motion is driven by other atoms in the system with which 
hyperfine interaction strengths 8' are larger than S with 
the spins of interest if r^ 1 ~ 6'. There seems to be some 



indirect experimental evidence that > Th, or, equiva- 
lent^, that ThSh < 1 in Fes and in Mni2 clusters. Resolu- 
tion of the tunneling pair of energy eigenstates has been 
reported in FesliS and in Mn^Jij under large transverse 
fields, both at 680 MHz. This implies that r > 1CT 9 
seconds in both cases. This is to be compared with the 
value of 1/Sh that would equal to [see, Eqs. (||) and 
©] if T h Sh > 1 were fulfilled. Now, Sh w 0.02K for 
Fes (as follows from measured hyperfine fields of 10 _3 T 
reported by WernsdorferE£l, g = 2, and S = 10), that 
is, USh w 0.5 x 10~ 9 s. For Mm 2 , Sh w 10~ 2 T, (see 
Ref.H), from which 1/Sh « 0.5 x 10~ 10 s follows. It fol- 
lows that, in fact, ThSh > 1 is not fulfilled, neither for 
Fes nor for Mni2, suggesting that tunneling in Fes and 
in Mni2 proceeds with some degree of coherence, even 
when no external transverse field is appliedJIS 
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